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Abstract. The Neutrinoless double beta Decay allows to determine the effectice Majorana 

electron neutrino mass. For this the following conditions have to be satisfied: 

(i) The neutrino must be a Majorana particle, i. e. identical to the antiparticle. 
i-^ . (ii) The half life has to be measured. 

(iii)The transition matrix element must be reliably calculated. 
Oj, (iv) The leading mechanism must be the light Majorana neutrino exchange. 

1^ ■ The present contribution studies the accuracy with which one can calculate by different 

^ ; methods:(l) Quasi-Particle Random Phase Approach (QRPA), (2) the Shell Model (SM), 

(3) the (before the variation) angular momentum projected Hartree-Fock-Bogoliubov method 

(PHFB)and the (4) Interacting Boson Approach (IB A). 
^ I In the second part we investigate how to determine experimentally the leading mechanism for 

QQ . the Neutrinoless Double Beta Decay. Is it (a) the light Majorana neutrino exchange as one 

^^ ' assumes to determine the effective Majorana neutrino mass, ist it the heavy left (b) or right 

V^ , handed (c) Majorana neutrino exchange allowed by left-right symmetric Grand Unified Theories 

CO ■ (gut's). Is it a mechanism due to Supersymmetry e.g. with gluino exchange and R-parity and 

»«f\ ' lepton number violating terms. 

("^ , At the end we assume, that Klapdor et al. [l] have indeed measured the Neutrinoless Double 

^SJ ' Beta Decay(, although heavily contested, )and that the light Majorana neutrino exchange is the 

leading mechanism. With our matrix elements we obtain then an effective Majorana neutrino 

mass of: 

k>J ■ < m„ >= 0.24[eV](exp ± 0.02; theor. ± 0.01) [el^] (1) 



1. Introduction 

The electron neutrino mass can be determined in the single beta decay (Mainz, Troisk and 
in the future KATRIN/Karlsruhe) and for Majorana neutrinos also in the neutrinoless double 
beta decay. Although with both methods one has not yet succeeded to determine the effective 
electron neutrino mass, the neutrinoless double beta decay has in the future the potential to 
to measure the effective Majorana neutrino mass even to smaller values than in the single beta 
decay. For this the neutrino needs to be a Majorana particle (identicall with the antiparticle), 
the transition matrix element must be able to be calculated reliably and the light Majorana 
neutrino exchange must be the leading mechanisme for the neutrinoless double beta decay. 

^ Talk at the Schladming meeting, February 25th. to March 3rd. 2012. 



The purpose of this contribution is first to investigate the accuracy, with which one can calcu- 
late the transition matrix elements. Second we want to find an experimental proceedure, which 
allows to determine the leading mechanisme for the neutrinoless decay: Is it the light Major ana 
neutrino exchange as usually assumed or the exchange of a left or right handed heavy neutrino 
allowed in left-right symmetric Grand Unified Theories (GUT's). An additional possible mech- 
anisme is due to suppersymmetry (SUSY) by a trilinear term in the Lagrangian, which violates 
lepton number and R-parity conservation. 



The transition matrix elements for the light Major ana neutrino exchange are presently 
calculated with four different methosd. They are reviewd in chapter, where they are compared 
with each other with their advantages and their drawbacks ^. The different methods are the 
Quasi-particle Random Phase Approximation (QRPA) for spherical (chapter 2) and deformed 
nuclei H S H El El [7] , the SheU Model (SM) El HO] , the Projected Hartree-Fock-Bogoliubov 
(PHFB) approach |1H [T2l [T3l [T^ [T5] (In reference jl43 a missing factor two in all older works 
of the group has been pointed out.) and the Interacting Boson Model (IBM2) |16j . 
In the third chapter we do not assume, that the light left handed Majorana neutrino exchange 
is the leading mechanism |17 1 I18 1 [T9]. We search for a possible experimental identification of the 
leading mechanism even in cases, where two or even three equally strong mechanisms interfere. 
The diagram for the neutrinoless double beta decay with the exchange of the light Majorana 
neutrino is shown in figured] for the decay of 32^644 to 34 £'642 through the intermediate nucleus 
33^843. 




Figure 1. The light Majorana neutrino exchange mechanism for the neutrinoless double beta 



decay of 32^644 to 345642 through the intermediate nucleus 33^^43 is shown, the neutrino 
must be identicaql with the antineutrino (Majorana neutrino) and helicity must not be a good 
quantum number, this means the neutrino must be a massive Majorana particle. 



In addition to the light left handed Majorana neutrino exchange, one has other possible 
mechanisms as cause for the neutrinoless double beta decay: Grand Unification (GUT), 
Supersymmerty (SUSY) and extensions to extra dimensions. We shall discuss here extension to 
GUT'S and SUSY. 
In a left-right symmetric model of GUT we have: 



Wi = cos '&GUT + sin '&GUTW2 = — sin '&GUT + cos "dcuT 



(2) 



Wi is the usual vector boson of Rubbia and coworkers of 80.4 GeV mainly responsible for 
the left handed weak interaction. The mixture in eq. [fallows at each vertex in figure [1] left and 
right handed interactions. In addition one has SUSY contributions, where mainly the threelinear 
terms are responsible for the lepton number violation [T71 [181 13 EO]- This will be discussed in 
chapter three of this contribution. 
Fermis Golden Rule of second order time dependent perturbation theory yields: 

-7-0^^ _ /" ^F c ^ \^ < f\Hw\k >< k\Hw\i > . s 



T^^ = Ml""- <mu> +M^ < tani9 > +Mwn < i^f > + 

M2 

MsusY ■ Afn + Mnr < t^ > +••• (4) 

The first term on the right hand side of eqn. (HD with the matrix element M^ and the 
effective Majorana mass 

<m^>= Y^ {Uekf • mfc^ = Y^ e^^^fe . \Uek? ■ rriky (5) 

fc=l,2,3 fc=l,2,3 

with 

^e= Y ^(^k^k (6) 

fc=l,2,3 

is the light Majorana neutrino exchange contribution. 

2. The different Many Body Approaches for the Oi^/3/3 Matrix Elements. 

The Quasiparticle Random Phase Approximation (QRPA) is used by the groups in Tiibingen, 
Bratislava and Jyvaskyla [3 El H El El Ej, while the Strasbourg-Madrid group [SI El [TO] 
uses the Shell Model (SM), Tomoda, Faessler, Schmid and Gruemmer [11] and Rath and 
coworkers [12l [T3l ^^ use the angular momentum projected Hartree-Fock-Bogoliubov method 
(HFB), and Barea and lachello [16] use the Interacting Boson Model (IBM2, which distinguishes 
between protons and neutrons). 

The QRPA [H [3l H] [5] has the advantage to allow to use a large single-particle basis. Thus, 
one is able to include to each single nucleon state in the QRPA model space also the spin-orbit 
partner, which guarantees that the Ikeda sum rule [21] is fulfilled. This is essential to describe 
correctly the Gamow- Teller strength. The SM [8] is presently still restricted to a nuclear basis 
of four to five single-particle levels for the description of the neutrinoless double beta decay. 
Therefore, not all spin-orbit partners can be included and, as a result, the Ikeda sum rule is 
violated by 34 to 50% depending on the single particle basis used. 

In QRPA one starts from the transformation to Bogoliubov quasiparticles : 

of = Uicl - ViCi- (7) 

The creation c- and annihilation operators of time reversed single-particle states q are usually 
defined with respect to oscillator wave functions [3j. The single-particle energies are calculated 
with a Woods Saxon potential ^. 



The excited states \m) with angular momentum J in the intermediate odd-odd mass nucleus 
are created from the correlated initial and final 0^ ground states by a proton-neutron phonon 
creation operator: 

\m) = Ql\0+)- HQl\0+)=EmQl\0^). (8) 

Ql = Y.i^^^l-^^^-], (9) 

a 

which is defined as a linear superposition of creation operators of proton-neutron quasiparticle 
pairs: 

Ai = [aWkijM, (10) 

For the present presentation the complication of angular momentum coupling, which must 
and is included in the quantitative calculations, is not shown. 

The inverse 0^'/3/3 lifetime for the light Majorana neutrino exchange mechanism is given as 
the product of three factors, 

T^J^y^ = G^'' iMl-'f ■< m^ >2 (11) 

where G^'^ is a calculable phase space factor, A^^J'^ is the Oiyl3/3 nuclear matrix element, and 
< mu > is the (nucleus-independent) "effective Majorana neutrino mass" ([5]). 

The expressions for the matrix elements M^'^ and the corresponding Ou/S/S transition 
operators are given, e.g., in Ref. ^: 

9A 

The SM approach has been applied by the Strasbourg-Madrid group |10] to neutrinoless 
double beta decay [8] using the closure relation with an averaged energy denominator. In this 
way one does not need to calculate the states in the odd-odd intermediate nuclei. The quality 
of the results depends then on the description of the O"*" ground states in the initial and final 
nuclei of the double beta decay system, e.g. Ge — )-^^Se, on the nucleon-nucleon interaction 
matrix elements fitted by the Oslo group in neighbouring muclei and on the average energy 
denominator chosen (fitted) for closure. The Oi//3/3 transition matrix element ()12p simplifies as 
shown in equations (5) to (11) of Ref. [9|. Since the number of many body configurations is 
increasing drastically with the single-particle basis, one is forced to restrict for mass numbers 
A = 7G and A = 82 in the SM the single-particle basis to lps/2, O/5/2, lpi/2 and 0^9/2- In ref. [9] 
the ^^Se nucleus is calculated in the SM for five basis single-particle levels including also O/7/2. 
For the mass region around A = 130 the SM basis is restricted to 0(77/2, 1^3/25 1^5/2, 2si/2 and 
0/111/2 levels. The problem with these small basis sets is that the spin-orbit partners O/7/2 and 
Og'7/2 have to be omitted [9j. The SM results then automatically violate the Ikeda Sum Rule 
(ISR) [2l], while the QRPA satisfies it exactly. The Ikeda sum rule is: 

S_-S+ = 3{N - Z), (13) 

A A 

M k I 

For 5+ the subscripts at the isospin rising and lowering operators are exchanged. 

Figure [2] shows the QRPA contributions of different angular momenta of the neutron pairs, 
which are changed in proton pairs with the same angular momenta. In figure [2] the left bar is 
the result for ^^Se obtained with the single-particle basis 1^3/2, O/5/2, lpi/2 and 0^9/2 used in the 
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Figure 2. (Color online) Neutron pairs with different angular momenta J'^ transform into 
proton pairs with the same angular momenta. The different contributions for the angular 
momenta to the total M^^ calculated within the QRPA and different basis sizes for the Ou^P 
decay ^^Se— >-^^Kr is shown. The left bar is calculated with the same basis of four levels, 
lP3/2i 0/5/2) lpi/2 and 0^9/2, used in the shell model calculations [HI O [10]. The Ikeda Sum 
Rule (ISR) [21] is exhausted by 50%. The second bar from the left includes in addition the 
1/7/2 level, one of the two missing spin-orbit partners given for the ®^Se nucleus in ref. [S] 
for the shell model. The ISR is exhausted by 66%. The third bar from the left includes 
both missing spin-orbit partners O/7/2 and O57/2 amounting in total to 6 single-particle levels. 
The ISR is fulfilled by 100%. This leads to the increase in the neutrinoless matrix element 
from 1.12 to 4.07. The right bar represents the QRPA result with 9 single-particle levels 
(1/7/2, 2p3/2, 1/5/2, 2pi/2, 1^9/2 > 2^5/2, 3si/2, 2^3/2, 157/2-)- The matrix element gets only slightly 
increased to 4.27. The spin-orbit partners are essential to fulfill the Ikeda Sum Rule (ISR). 



SM. The ISR is exhausted by 50% only. The second bar from the left represents the result with 
addition of the I/7/2 level. The ISR is exhausted by 66%. The third bar from the left shows the 
result obtained by inclusion of both spin-orbit partners O/7/2 and Og'7/2 missing in the four level 
basis of the SM. The ISR is 100% fulfilled. For the right bar the basis is increased to 9 single- 
particle levels for neutrons and protons (O/7/2, IP3/2, O/5/2, lpi/2, 0^9/2, 1^5/2, 2si/2, 1^3/2, O57/2). 



In the last ten years P. K. Rath and coworkers |12^ \T3\ [T^ have pubUshed a whole series 
of papers (see references in Ref. [13J) on 2vj3l3 decay and, since 2008, also on {)vj3j3 decay, in 
which they used a simple pairing plus quadrupole many body Hamiltonian of the Kumar and 
Baranger type [22j to calculate the neutrinoless double beta decay transition matrix elements 
with angular momentum projection from a Hartree-Fock-Bogoliubov (HFB) wave function after 
variation. Schmid [23j did show, that with the assumption of a real Bogoliubov transformation 
(real coefficient A and B), axial symmetry 

M 

ai = ^{Aiac\ + Bi^a) (15) 

and no parity mixing, only 0^, 2"'",4"'", . . . . nucleon pairs and excited states are allowed (See 
eqn. (4.2.3) on page 603 of ref. [23]). Rodriguez and Martinez-Pinedo start with the projected 
HFB approach but allow admixtures of different deformations using the Generator Coordinate 
Method (GCM) and the Gogny force [23]. The QRPA and the SM do not have this restriction 
like PHFB and also in the PHFB with the deformation GCM extension (GCM-PNAMP) [15] . 
Figure [3] shows on the other side, that the projected HFB approach is restricted to con- 
tributions of neutron pairs with angular momenta 0"*" , 2"'" , 4"*" , . . . . In addition, one sees that 
the contributions of transition of higher angular momentum neutron to proton pairs 2"'" , 4"*" , . . . 
are drastically reduced compared to the QRPA and the SM see fig. [3l The reason for this is 
obvious: in a spherical nucleus the HFB solution contains only seniority zero and no stronger 
higher angular momentum pairs. The double beta decay system 32^644 — >-34 ^642 has only small 
deformations and thus a projected HFB state is not able to describe an appreciable admixture of 
higher angular momentum pairs for 0"*" — >■ 0"*" transitions as can be seen in ref. [23] . The higher 
angular momentum contributions increase drastically with increasing intrinsic quadrupole and 
hexadecapole deformations of the HFB solution. 

The results in figure [3] are calculated in.Tuebingen by K.W. Schmid [23] within the HFB with 
angular momentum and particle number projection before variation with an improved Gogny 
force p4J adjusted in a global fit to properties of many nuclei. 

To have also l'*',3~'', 5"*", . . . neutron pairs contributing one has to use a Bogoliubov 
transformation with complex coefficients A and B ()15p . To have also 0~, 1~, 2~, 3~, 4~, 5~, . . . 
one has to allow parity mixing in the Bogoliubov transformation (|15p. But even allowing 
all different types of angular momentum and parity pairs one would still have an unnatural 
suppression of the higher angular momenta especially for smaller deformations. This handicap 
could probably be overcome by a multi-configuration HFB wave function [23j with complex 
coefficients and parity mixing in the Bogoliubov (|15p transformation. 

The IBM (Interacting Boson Model) [TB] can only change 0"*" (S) and 2^ (D) fermionic pairs 
from two neutrons into two protons. In the bosonization to higher orders this leads to the 
creation and annihilation of up to three "s" and "d" boson annihilation and creation operators 
in Ref. [TB]. But all these terms of equation (18) of reference [16] originate from the annihilation 
of a 0"*" (S) or a 2"'" (D) neutron pair into a corresponding proton pair with the same angular 
momentum. The higher boson terms try only to fulfill the Fermi commutation relations of the 
original nucleon pairs up to third order. The IBM can therefore change only a 0"*" or a 2"'" 
neutron pair into a corresponding proton pair. 

We have also calculated [6]:7] the transition matrix elements of the light left handed Majorana 
neutrino exchange M^'^ including deformations. Different deformations for the inital and the 
final nuclei are allowed. The BGS overlaps are included. The quadrupole deformations are taken 
from the reorientation Coulomb excitation of the 2"'" states. 
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Figure 3. Neutron pairs with different angular momenta contributing to tfie neutrinoless double 
beta decay transition matrix elements for '^^Ge ^'^^Se calculated from a HFB wave function 
with angular momentum and proton and neutron particle number projection before variation. 
The Fermi, the Gamow- Teller and the total contribution including the tensor part as defined 
in eq. ()12p are separately given. The nucleon-nucleon interaction is an improved Gogny type 
force [i24J. The deformations /Sce = —0.08 and I3se = 0.11 correspond to the minima of the 
projected HFB total energy. The results for the transition matrix elements are qualitatively and 
almost quantitatively the same for the experimental deformation from the Coulomb reorientation 
effect: /Sce = 0.16, /35e = 0.10 and also for different forces. The angular momenta of the neutron 
pairs are in the PHFB approach with axial symmetry, real coefficients and no parity mixing 
restricted to 0^ , 2^ , 4"*" , . . . . In addition the contributions of higher angular momentum neutron 
pairs 2"*", 4+, . . . are drastically reduced compared to the QRPA and the SM. 
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-^charge 

The deformation reduces the matrix elements [61I3 in Ge slightly from 5.30 to 4.69 by 10% 
only. This is within the error of the matrix elements (see figure [5]) . But the reduction of the 
matrix elements is severe in strongly deformed systems with different deformations for the initial 
and the final nuclei. In the system ^^^Nd -^^^^ Sm the matrix element is reduced from 6.12 to 



Table 1. Nuclear matrix elements (NME) for deforemed "1" and for spherical "0" QRPA 
calculations Ml"" for Oi//3/3 decays ^'^Ge^'^^Se, i5°Nd^i5°Sm, 160q^^160 j-,y r^^ie BCS overlaps 
are taken into account. In the last two columns the Ou/]/] matrix element A^J^" and the half-lives 
for assumed < rui, >= 50 meV are shown. 



A 


Def. 


gA 


M^-^ 


^1% • [W^'y] 
< rriy >=50 meV 


76 


it-l ■>■> 

"0" 


1.25 
1.25 


4.69 
5.30 


7.15 
5.60 


150 


"1" 
"0" 


1.25 
1.25 


3.34 
6.12 


0.41 
0.12 


160 


"1" 


1.25 


3.76 


2.26 



3.34 (see table 2) and in the strongly deformed system ^^^Gd — s-i^o j^y ^ggg ^^ble 1) one obtains 
a matrix element of 3.76 (see table 2). The single nucleon basis in these deformed calculations 
are determined in a deformed Woods-Saxon potential. The results are then expanded into a 
deformed oscillator basis with the same deformation parameter and the appropriate oscillator 
length in seven oscillator shells [U [7]. The deformed result for "^^^Nd — ^^^^ Sm is included in 
figure O 

3. How to find the Leading Mechanisms for the Neutrinoless Double Beta Decay? 

Normally one assumes, that the first term of eq. (JH) is the leading one and with the experimental 
data and the matrix element for the light left handed Majorana neutrino exchange M.^ one 
can determine the effective Majorana neutrino mass (0). But in Grand Unification (GUT) and 
Supersymmetry (SUSY) ad^'^'^^""' mQr.i^oTiioT,io fr^v +i^q n^nifT-^vir^iogg Double Beta Decay (Quj3j3) 
are possible. 



W, 



W, 



Figure 4. Light Majorana neutrino exchange. The two neutrons in the initial nucleus in 
the ground state 0"*", which change into two protons are are coming from the left and are 
characterized by the two down quarks d in the two neutrons, which change into two up quarks 
u in two protons in the final nucleus. The vector bosons Wl mediating the left handed weak 
interaction are coupled by the expansion coefficients C/efe ( ED to the neutrino mass eigenstates 
rriky 



The matrix element of this gold plated term is proportional to: 

M'iriUght i^l) ^ Yl K,-PL-^^PL-Kk = \ E e^-^-lf^r^P^ifc. (17) 

fc=l,2,3 ^ ^^'^ '^ A;=l,2,3 

The exchange of a heavy left handed Major ana neutrino: 



N,= Y. ^ekNk^ Y. ^"''\Ul\-Nu (18) 

fc=l,...6 fc=4,5,6 

with a^ the Majorana phases for these heavy left handed Majorana neutrinos. 
M'^'^iheavyNr.)^ Y ^ek ' Pl , _\ Pl ■ U^k = - Y e'^''^ ' W^kl'/MkN (19) 

fc=4,5,6 ^ '^^ fc=4,5,6 

The lepton number and R-parity violating contributions in SUSY are the trilinear terms and 
the coupling of the lepton superfields to the Higgs particle. 

Wi^ = hjk ■ Li ■ Lj ■ El + A'yfc • Li ■ Qj ■Dl + fii-U-H2 (20) 

The lepton L, E and quark Q, D left (L) and right (R) handed superfields are defined as: 



kL ^^^ kL 

The indices i,j,k run over the three families for leptons: e,/i, r and for quarks: d,s,b. The 
subscripts L and R characterize left and right handed superfields. The tilde indicates SUSY 
particles like selectrons , sneutrinos and squarks. 
The inverse half life is given by: 



T^J^ ln2 "^ '^'''^' '^■'"'"'^ 
tinlM'^nl + rix'Mlff + \vnr\'\M'^r\'\ (22) 

with: 

^^ = <'^-> = (1^- |2 . ^^ ^ ^2^a,^ . |f^.^|2 . ^^ ^ ^2^as, . ^jju^^2 . ^^)/^^ (23) 

7np_ 



riATT - lr/^|2 . "^P -L p2i«54 . |rrAf|2 . "^P , 2iae4 . Ir/A^p . ^P (oA) 

riNL-\U,^\ ^^^+e \U,,\ ^^^+e \U,,\ ^^^ (24) 



"ifc = Oi - ttfc (25) 



We restrict here for SUSY to the trilinear terms ()20p .Thev are lepton number and R parity 
violating. They can contribute by gluino or by neutralino exchange. 
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Figure 5. (Color online) Transition matrix elements of the neutrinoless double beta decay for 
the different approaches: QRPA [3111], the SM (HJ [IHl E] , the projected HFB method [H], the 
projected HFB with the Generator Coordinate Method (GCM) with deformations [15] (GCM- 
PNAMP) and the IBM2 [IS]. The error bars of the filled circles for the QRPA are calculated 
as the highest and the lowest values for three different single-particle basis sets, two forces 
(Bonn CD and Argonne V18) two different axial charges gA = 1-25 and the quenched value 
gA = 1-00 and two different treatments of short range correlations (Jastrow-like [2S] and the 
Unitary Correlator Operator Method (UCOM) [2B]). The radius parameter is as in this whole 
work ro = 1.2 fm. The triangle with the tip up are the SM results [51 [TUl [U]. The triangle with 
the tip down represent the transition matrix element of the Interacting Boson Model 2 (IBM2) 
[T6] . The squares have been calculated by Pradfulla Rath and coworkers [H] with the correction 
of the factor 2 from December 2010 included, with which all previous results of Rath et al. have 
to be multiphed [121 [131 [Hj. The star (GCM-PNAMP) is a projected HFB calculation with 
the Gogny force [23] by Rodriguez and Martinez-Pinedo [15] allowing for different deformations 
with the Generator Coordinate Method (GCM). 



To test, if the light left handed Majorana neutrino exchange is the leading mechanism, one meeds 
at least experimental data of the neutrinoless double beta decay and reliable transition matrix 
elements in two systems. The light Majorana neutrino exchange is represented by the gold plated 
term, which is the first on the right hand side of eq.©. If the measurements and the matrix 
element are reliable enough [TU] and the light Majorana neutrino exchange is indeed the leading 
mechanism, one should in both and also in all other systems obtain the same effective Majorana 



neutrino mass. If two mechanisms are at the same time contributing, one must distinguish 
between non-interfering and between interfering mechanisms. The hght left handed Majorana 
neutrino and the heavy right handed neutrino exchange have neghgible interference |18j . 



Tf;if^-G^^{E,,Z) 



riu\'{MZY + \mR\\M%nY (26) 



To determine the absolute values of the two strength parameters rj^ and r/^vi? one needs at 
least two decay systems i. To verify, that these are indeed the leading mechanisms one needs at 
least a measurement in one additional decay system i. But if one forms ratios of half lives using 
the Tuebingen matrix elements for the light Majorana neutrino exchange and the heavy right 
handed neutrino exchange one otaines a very restricted allowed interval for these ratios [18] . 

0.15 < ; ,^^ -<0.18; 0.17 < ' „^ r < 0.22; 

1.14 < '\,^^ < 1.24: (27) 

The dependence of this ratios on the different parameters of the nuclear structure calculation 
is very minor. Due to the ratios the dependence on most changes drop approximately out. The 
ratios (|27p are calculated for the axial charge qa = 1-25. But the quenching of this value to 
QA = 1-00 has only a minor effect |18j . 

If the two leading mechanisms like the light Majorana neutrino exchange and the SUSY 
mechanism with gluino or neutralino exchange can interfer, the situation is a bit more 
complicated: Let us assume the relative phase angle of the complex strength parameters rj^ 
and ri\i is "Sy^x' . The inverse of the half life time the phase space factor G^,2{Eq, Z) is then: 



cos??,,v|r?.|-|??A'|--A^°;--^°A'; (28) 

One needs three decay systems to determine the absolute values of the parameters r/jy, r]x' and 
the relative phase angle i^iy^x'. At least one additional system is needed to verify, that indeed 
these two mechanisms are the leading ones. Again the ratios of the half lives are allowed to lie 
only in narrow regions [18]. If this is not the case, the chosen mechanisms are not the leading 
ones [18] . With CP conservation the strength parameters r] must be real and thus the relative 
phase angle is zero or 180 degrees. So the determination of '&u,x' allows to test CP conservation 
or violation. 

4. The effective Majorana Neutrino Mass. 

Before we summarize the results let us assume Klapdor-Kleingrothaus et al. p] have indeed 
measured the neutrinoless double beta decay in "^^Ge, although the general belief is, that this 
still needs confirmation. From the half life given by Klapdor et al. [1] one can derive with our 
matrix elements the effective Majorana neutrino mass ([5]). 

T^J^C'^Ge, exp Klapdor) = (2.23 + 0.44 - 0.31) • 10^^ [years]; (29) 



With our matrix elements one obtains the effective hght left handed Majorana neutrino mass 
under the assumption, that the light Majorana exchange is the leading mechanism. 



< niy >= 0.24[ey] (expzb 0.02; t/ieor. ±0.01) [eF] (30) 

The uncertainty (error) from experiment is 0.02 [eV], while the theoretical error originates 
from the uncertainties of the QRPA matrix elements as indicated in figure [5j The theoretical 
error is 0.01 [eV]. 

5. Conclusions 

Let us now summarize the results of this contribution: 

The Sheh Model (SM) [8l|10l|9] is in principle the best method to calculate the nuclear matrix 
elements for the neutrinoless double beta decay. But due to the restricted single-particle basis 
it has a severe handicap. The matrix elements in the ^^Ge region are by a factor 2 smaller than 
the results of the Quasiparticle Random Phase Approximation (QRPA) [2l[3llH[5], the projected 
Hartree Fock Bogoliubov approach [14, 15J and the Interacting Boson Model (IBM2) [TB]. With 
the same restricted basis as used by the SM the QRPA obtains roughly the same results as the 
SM (figure [2]) , but the Ikeda sum rule [21j gets strongly violated due to the missing spin-orbit 
partners in the SM single-particle basis. 

The angular momentum projected Hartee-Fock-Bogoliubov (HFB) method [T^ is restricted in 
its scope. With a real Bogoliubov transformation without parity mixing and with axial symmetry 
(jlSp one can only describe neutron pairs with angular momenta and parity O"*", 2"*",4+,6^, . . . 
changing into two protons for ground state-to-ground state transitions. The restriction for the 
Interacting Boson Model (IBM) |16] is even more severe: one is restricted to 0"^ and 2"'" neutron 
pairs changing into two protons. 

A comparison of the Oz^/3/3 transition matrix elements calculated recently in the different many 
body methods: QRPA with realistic forces (CD Bonn, Argonne V18), SM with nucleon-nucleon 
matrix elements fitted in neughbouring nuclei, projected HFB [14j with pairing plus quadrupole 
force f22], projected HFB with the deformation as Generator Coordinate (GCM+PNAMP) |15j 
and with the Gogny force[24] and IBM2 [16] is shown in Fig. [5j 
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